A Metal-Insulator transition induced by Random Dipoles 
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We study the localization properties of a test dipole feeling the disordered potential induced by 
dipolar impurities trapped at random positions in an optical lattice. This random potential is marked 
by correlations which are a convolution of short-range and long-range ones. We show that when 
short-range correlations are dominant, extended states can appear in the spectrum. The effect of 
long-range correlations is to restore localization over the whole spectrum and to enhance or reduce 
localization at the band center, depending on the localization regime determined by short-range 
correlations. 
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I. INTRODUCTION 

Interference effects induced by random potentials 
deeply modify wave diffusion unto even stop wave trans- 
port on a length Cioc^ the localization length. This phe- 
nomenon called strong localization or Anderson localiza- 
tion [l|, can be observed in classical waves such as 
acoustic waves microwaves light Q, as well as 
in quantum waves, such as electronic [611 or matter waves 
in real or in momentum space [111 [T2j . Anderson 

theory relies essentially on two main assumptions: (i) the 
potential has to be (5-correlated; (ii) the wave has to be 
noninteracting. If these conditions are fulfilled, Anderson 
localization occurs always in one (ID) and two dimen- 
sions (2D), and depending on the disorder strength and 
the energy in three dimensions (3D) [l^ . A localization- 
delocalization transition can appear in low dimensions 
if correlations and/or interactions are introduced in the 
systems. If on the one hand noninteracting waves exist 
in nature, on the other hand real uncorrelated poten- 
tials do not exist. Correlated disordered potentials can 
be roughly divided into three wide classes, accordingly to 
the behaviour of the two-point auto-correlation function 
C(£) and of its Fourier transform S{k). The first class is 
marked by a correlation function C{i) which decays expo- 
nentially on a length £ with 'k /kmax < £ <^ L, L being the 
length of the system and kmax the largest wavevector al- 
lowed by the system. These potentials, called short-range 
correlated potentials, may introduce resonance energies 
in the spectrum inducing delocalization of a significant 
subset of the eigenstates. This happens for instance in 
the random-dimer model (RDM) and in its dual counter- 
part (DRDM) [ll-[l3, 

in which the sites of a lattice are 
assigned energies Sa or e?, at random, with the additional 
constraint that sites of energy Sb always appear in pairs 
(RDM) or never appear as neighbors (DRDM). The sec- 
ond class of correlated potentials is marked by a spectral 
function S{k) that is non-zero in a finite fc-range. This 
is, for instance, the case of the speakle 18, 19] where 
C{£) is a sine function. For this kind of potentials, there 
exists a critical energy at which the localization length 




FIG. 1: (Color online) Schematic representation of the phys- 
ical model. Dipolar impurities (green spheres) are trapped at 
the minima of the optical lattice, occupying random positions. 
The test dipole (red sphere), excited to another internal level, 
feels a shallower optical potential and a disordered potential 
due to the dipolar interaction with the impurities. 



increases abruptly mimicking the presence of a mobility 
edge in finite-size systems [20, [21|. In the third class 
C{£) decreases algebraically as ~ 1/^'^, and both C{£) 
and S{k) are nonzero over the whole real and k spaces. 
In practice there are no length scales characterizing the 
disorder, which is scale-free. These potentials are com- 
monly called long-range correlated. In this case, it has 
been observed that correlations can have different effects 
depending on the energy region under consideration. In 
particular, for discrete models, a reduction of the local- 
ization length has been observed at the band edges and 
conversely an enhancement at the band center has been 
reported [l^]. In this context also the presence of mo- 
bility edges has been claimed [231 . although these results 
stirred some controversy [l^, l25j |. 

Very often, especially in the case of long-range correla- 
tions, these studies rely on toy models characterized by 
ad-hoc correlation functions, creating almost no connec- 
tion with possible experimental implementations. In this 
paper we propose a physical model for a random potential 
where long-range and short-range correlations arise nat- 
urally from the system itself, exploiting the properties of 
the dipolar interaction [26,, J7,] in ultracold atomic (28ri3l| 
or molecular [H, [11] gases. The model considers a series 
of dipoles pinned at random positions at the minima of 
a deep optical lattice. The dipoles are polarized perpen- 
dicularly to the lattice axis, so that dipole-dipole inter- 
action is repulsive. In this way, for low enough densities, 
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there are no double occupancies, and, if the dipole-dipole 
interaction is large enough, consecutive occupations are 
avoided too. This set of trapped dipoles, referred as im- 
purities, will create a disordered potential [l^ [sj, HI] 
for another dipole, the test dipole, which is excited to a 
different internal level and is able to move through the 
lattice (see Fig. [1]). Short-range correlations arise from 
the distribution of the impurities, while long-range cor- 
relations are due to the dipolar interaction between the 
test dipole and the impurities. 

We study the localization properties of the test dipole 
in the correlated potential realized by the impurities, 
highlighting the role played by short and long-range cor- 
relations. In particular, depending on the parameters of 
the model, we observe that short-range correlations can 
introduce a discrete set of extended states in the system 
while long-range correlations tend to restore localization 
and lead to counterintuitive effects on the localization 
length of the system. 

The paper is organized as follows. The model is pre- 
sented in detail in Sect. HI] and the Hamiltonian for the 
test dipole is derived. In Sect. IIIH we study the localiza- 
tion properties of the model by using a renormalization- 
decimation scheme for the calculation of the localization 
length. Then, a detailed discussion of the role played 
by short and long correlations is presented in Sect. IIVI 
Finally, we draw our conclusions in Sect. [V] 

II. THE PHYSICAL MODEL 

We consider a very dilute gas of dipolar impurities 
trapped in a deep ID optical lattice and a single test 
dipole excited to a different internal level with the same 
dipole moment D. We superimpose to the optical lattice 
a very elongated cigar-shaped harmonic confinement, so 
that we can practically neglect the weak axial harmonic 
confinement along the z axis. We assume both for the 
impurities and the test dipoles frozen radial dynamics 
into the lowest state of the radial harmonic confinement 

0c.(rx) = ^e-'-i/^'^S with a. = 

(1) 



where a; = o;^^-* and uj = uj^'^^ for impurities and test 
particle respectively. 

Furthermore, assuming that the optical lattice felt by 
the impurities is so strong to completely freeze their dy- 
namics along z, we can describe the motion of the test 
particle of mass m along the lattice axis z by the follow- 
ing ID effective Hamiltonian 



H = -^-^ + sEn^sm''ikz) + Vd{z), (2) 

where Er — ir^h^ /2md\ is the recoil energy, dh the lat- 
tice spacing of an optical lattice generated by a laser of 
wavelength = 2dL, and s the adimensional lattice 
depth. 

The random potential Vd{z) results from the dipolar 
interaction of the test particle with the impurities pinned 
in the lattice. It is given by the convolution of the density 
distribution of the impurities p{z) (see Eq. [4|) and the 
effective one dimensional dipolar interaction U\^{z) (see 
Eq.©: 



Vd{z)^ Jdz'piz'M^iz-z'). (3) 

The density distribution p(z) is given by the sum of the 
Wannier functions w^-^\z) localized around the sites I 
occupied by the impurities 

p(z) = X^M^)(z-rdi)l^ (4) 

The random distribution of the occupied sites I intro- 
duces disorder in the system. 

The effective one dimensional dipolar potential U^^{z) 
is obtained after integration of the dipolar interaction 
Udd{r) = D^[l - 3{D ■ r)^]/|f|^ in the radial directions 
and is given by [s^ 



mdi^)^ / d^'±rfr_,|0^u,(r'i)n0^<T,(ri)|^{/dd(f-r') 
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where a± = J{crlii-, + (t'I^t))/'^, ol is the angle between the dipole D and the z axis and Erfc(z) is the comple- 
* mentary error function 
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Erfc(z) 



dt. 



(6) 



Note that the final expression that we obtained for 
U^^(z) is composed by two parts: a Dirac delta term 
at z = and a slowly decaying part. One can show 
that at large distances \z\ 3> ct±, the slowly decaying 
part reproduces the typical behaviour of the dipolar in- 
teraction, namely a decrease with the inverse cubic dis- 
tance U]2{z) ~ £)^(1 — Scos^ a)/|zp. For simplicity, in 
the present paper, we do not include contact interactions 
with the underlying idea that they can be switched off 
by exploiting Feshbach resonances [131 ■ However, they 
would simply modify the strength of the 5 part of the 
ID potential and provide an additional way to tune the 
parameters of the system. 

The tight binding form of the Hamiltonian is obtained 
using as a basis the set of Wannier states Wn{z) for the 
test particle. For the case of a single impurity pinned at 
site / we obtain 
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FIG. 2: (Color online) Site energies Aq, Ai and A2, and hop- 
ping energy 9 as function of a± in units oi cLl, for the case 
with 



of Dysprosium atoms with /d\ — 0.016i5ji, a = 7r/2, 
Al = 2dL ~ 570 nm, s^t) = 6 and = 30. The col- 
ored vertical lines labeled by different letters represent the 
set of values that we used for the calculation of the localiza- 
tion properties of the system. 



Hi = '^-J (I Wn){Wn+l I + I Wn+l){Wn |) 

71 

-J'^{\wi}{wi±i\ + \wi±i){wi\) (7) 

I Wn){Wn I . 

In Eq. ([7]), beyond the standard nearest neighbour tun- 
neling term J, we have included two terms due to the 
dipolar interaction: the first represents a nearest neigh- 
bour dipolar assisted hopping J'', while the second con- 
tains the on-site energies w^'l; at site n. The Hamiltonian 
parameters can be calculated using the following expres- 
sions 



We note that, for this choice of parameters, we can 
reasonably set 6* = 1 {J'^ — 0) and approximate the onsite 
energies for — ^| > 2 by the asymptotic expression 
A„_/ = X/\n-l\^ with A ^ D^/{Jdl). They are therefore 
independent on the value of a±. Also Ai does not depend 
significantly on the radial confinement. Conversely the 
value of Ao strongly depends on a± and can even vanish 
and become negative due to the anisotropy of the dipolar 
interaction. 

In the presence of several impurities the different con- 
tributions have to be included summing over the occu- 
pied sites I. Because of the impurity-impurity dipolar 
repulsions, we will impose that each impurity has to be 
preceded and followed by at least two empty sites. The 
resulting Hamiltonian, 
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,dd 



= I \wUz)\'\wP{z')\'U',^{z-z')dzdz'. 



The function u'^_i simply depends on the distance |h — Z| 
between the test particle and the impurity and provides 
the dipolar interaction between a single impurity and the 
test particle in the discretized formalism. 

In Fig. [21 we show the behaviour of the quantities 



= (J + J'')/ J and Xn-i -= 



,dd 



_i/J for \n-l\ = 0, f 
and 2, for D'^/dl = 0M6Er, a = tt/2, S(t) = 6 and 
S(/) = 30. This value of D"^ /d\ corresponds to the case 
of Dysprosium atoms trapped in an optical lattice gener- 
ated by a laser of wavelength = 570 nm [soj . 



H = - Wn){Wn+l I + I W„+l)(w„ |) 

n 

-f ^£„ I Wn){Wn |, 
n 

has hopping energies equal to J, and site energies 



(9) 



(10) 



The discretized impurity density distribution pi is a 
stochastic variable characterized by an average {pi) — C, 
where C is the impurity concentration, and by a correla- 
tion function Cp{i) — {pipi+i). The average and the cor- 
relation functions of the full potential can be extracted 
from the statistical properties of pi and the shape of the 
interaction potential. In particular, one can prove two 
important relations. First that the average value of the 
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full potential is simply given by the product of the im- 
purity concentration and the average value of the inter- 
action potential {udd) ■ 



,dd\ 



(11) 



second that 



namely that the two-point correlation function is given 
by the convolution of Cp{£) and 



(13) 



For the case of a dipolar potential lim^^oo Cu{i) oc £~'^, 
and for random impurities whose minimum distance is 
fixed to be three sites, as in our model, one has 



c 



1 - 2C 



e/2 



[Acos{K£) + Bsm{K£)] , (14) 



with K = ^ 



dl'[n - atan^(4-9C)/C], A = C - 



and 



2C) + {C- C^) cos(KdL)]/ sin(KdL). 
Thus we can conclude that the impurity distribution 
introduces short-range correlations, while the shape of 
the interaction u"^ is responsible for long-range correla- 
tions. The role and the competition between these two 
effects will be extensively discussed in Sect. IIVI 

For the full potential, let us also introduce the reduced 
correlation function, defined as 



Ce{£) 



(15) 



and the associated spectral density 
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Gn,n{E) 



Gi,n(E) 
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where Cioc{E) is the localization length, 
G{E) = {E — H)~^ is the Green's function related to the 
Hamiltonian H at energy E, and Gij{E) = {i\G{E)\j). 
The matrix elements Gi,n{E) and Gn,n, where N is 
the total number of lattice sites, have been computed 
by exploiting a renormalization-decimation scheme (38| . 
Our results are shown in the upper row of Fig. |3l where 
we consider increasing values of a±, corresponding to 
the vertical lines in Fig. [51 We use the same color code 
in all figures and label corresponding simulations with 
the same letters (a), (6), (c) and (d). Here and in the 
following, we fix C = 1/4, average over 100 configurations 
and consider system sizes up to 10^. As the value of 
a± is increased, we observe very different localization 
regimes. Notably for certain values of a±, we observe 
divergences of the localization length, corresponding 
to the appearance of metallic states in the spectrum. 
This suggests the presence of delocalization effects 
induced by the correlations of the physical model under 
consideration. 

More precisely, for large positive values of Aq, all states 
are clearly localized since the localization length is al- 
ways finite (first panel, (a)). By increasing a±, for van- 
ishing values of Aq the localization length exhibits two 
well defined peaks in two regions of the spectrum (sec- 
ond panel, (6)). Increasing a± further, corresponding to 
negative values of Aq, we observe the disappearance first 
of one of the two divergences (third panel, (c)), and then 
of both of them (fourth panel, (d)). In this last panel 
no divergences of the localization length are observed, 
but there are still peaks at the band edges, that recall 
the diverging behaviour previously observed. The dot- 
ted black lines correspond to the localization length Cl^^ 
evaluated in Born approximation, which corresponds to 
a second order perturbative calculation in the disorder 
strength iH-SlI, 



5(fc)=^ca^)e 



ike 



(16) 



In the following, we will use the square root of the 
variance of the full potential to quantify the potential 
strength W = ^(4) - i^n)^- 



III. NATURE OF THE SPECTRUM 

We study the nature of the spectrum of the tcst-dipole 
by evaluating the Lyapunov exponent A{E) through the 
asymptotic relation 



a^^\e) 



1 _ S{2k{E)) 



(18) 



where the connection between k{E) and the energy is 
given by the following relation E — (£„) -I- 2Jcos(fc(ii). 
Let us note that the Born approximation gives, by con- 
struction, a symmetric localization length around the av- 
erage value of the disorder (e„), since the spectral density 
S{k) associated to c{£) is always a symmetric function of 
k. Despite of this fact, there is a noticeable agreement 
between the Born approximation and the exact numerical 
results. 
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FIG. 3: (Color online) Top panel: localization length Cioc in units of the lattice spacing as a function of the energy in 
units of J for C — 1/4, system sizes up to 10^ and averaged over 100 configurations. The black dotted lines correspond to the 

(2) 

localization length calculated in Born approximation ClJ^. Bottom panel: reflection coefficient TZ of the single impurity as a 
function of the energy in units of J. The vertical dashed lines indicate the energies for which the reflection coefficient vanishes 
TZ{E) — 0. From left to right, different plots refer to increasing values of a±, corresponding to the vertical lines in Fig. [21 



IV. ROLE OF CORRELATIONS 

With the aim of understanding why we are observing 
the appearance and the disappearance of metaUic states 
in the spectrum by varying the radial confinement (and 
thus the ID dipolar interaction), we analyse the effects of 
the short-range correlations introduced by the impurity 
density distribution and of the long-range correlations in- 
troduced by the dipolar potential between the test dipole 
and the dipolar impurities. 

A. Short-range correlations 

In order to isolate the role of short-range correlations 
in our model, we calculate the reflection coefhcient for 
the case of a single impurity, as in Eq. ([7]), trapped at 
the site of index of an infinite lattice. To get analyti- 
cal results, we first neglect on-site contributions beyond 



nearest neighbours. Therefore we assume that a single 
dipolar impurity modifies just a trimer of on-site energies 
{Ai, Ao, Ai} and we generally assume that it can modify 
also the hopping with nearest neighbouring sites 6. 

A plane wave (n|fc) ~ ^ikud^ ^ eigen- 
state of the unperturbed periodic Hamiltonian 

= Y.t=-oo~Ji\ Wn){Wn+l \ + | W„+l)(w„ |) with 

energy E = — 2 Jcos(fcc?L), is perturbed by the impurity 
and results in the following wavefunction 



where r is the transmittance and r the reflectance. 
Using the scattering formalism combined to the 
renormalization-decimation scheme [l^, H^, H^, we ob- 
tain the following analytical formula for the reflection 
coefficient TZ = |rp of the single dipolar impurity 



n 



{^1 (f )' - f [1 - + ^? + -^i^o] - 2(?2Ai + Ao + A^Ao}' 



{^1 (f)' - f [1 - + ^? + ^i^o] - 202Ai + Ao + A^Ao}^ 
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(20) 
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In the bottom row of Fig. |31 we plot TZ for the same 
parameters used for the calculation of the localization 
length, i.e. 6 — 1 and Ao and Ai taken from the curves 
in Fig. [21 We observe that the calculation of the re- 
flection coefficient of the single impurity provides a very 
good understanding of the behaviour of the localization 
length: the energies where TZ tends to zero are exactly 
those where the localization length exhibits very large 
anomalous values. There is therefore a direct connec- 
tion between the appearance of metallic states in the 
spectrum and the scattering properties of the single im- 
purity. It has been previously shown by Dunlap et al. 
[14| that this kind of single impurity analysis can be 
used to interpret the transport properties of a system 
of size N where there are several randomly placed impu- 
rities. More precisely they proved that in such systems 
the number of single-particle states that show a metal- 
lic behaviour, being extended over the full system, is of 
the order of ^/N . Notably this number of delocalized 
states is large enough to induce transport in the system 
and initially localized wavepackets show a superdiffusive 
spreading in the disordered potential. This means that 
this type of extended states are detectable in typical ex- 
pansion experiments which can be performed with ultra- 
cold atomic gases 0, II] . 

It is remarkable that making use of the simple analyt- 
ical expression ()20|) we can predict the localization prop- 
erties of a rather complex system and the occurrence of 
metallic states in the spectrum. Studying the solutions 
of the equation 

TliE) = (21) 

as a function of Aq and Ai, one can extract the phase di- 
agram in Fig. m We identify four different regions in the 
diagram depending on the number of solutions of Eq. (|21|) 
and on their values. More precisely, if the solutions are 
both imaginary, no divergences are expected and all the 
states are exponentially localized (red region). If the so- 
lutions are real and inside the single impurity spectrum 
E = — 2 Jcos(fc(ii), namely when the roots satisfy the 
additional condition \E\ < 2 J, divergences are expected. 
Therefore, when the solutions are real, we can identify 
three additional scenarios: both solutions lie outside the 
spectrum (yellow region), one lies inside and one outside 
the spectrum (green region), both lie inside the spectrum 
(blue region). 

In the diagram shown in Fig. 21 we identify with mark- 
ers the values of Aq and Ai corresponding to the different 
plots of Fig. 131 where the dashed vertical lines mark the 
energies which verify the condition TZ{E) = 0. In partic- 
ular the square in the red region (a) corresponds to the 
first plot in Fig. |3] where all states are localized; the cir- 
cle in the blue region (6) corresponds to the second plot 
in Fig. |3l where we observe two resonances; the triangle 
in the green region (c) corresponds to the third plot in 
Fig.[3l where we observe one resonance; and the diamond 
in the yellow region (d) corresponds to the last plot in 
Fig. 131 where there are no resonances but the peak on the 
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FIG. 4: (Color online) Phase diagram induced by short- 
range correlations extracted from the reflection coeflicient in 
Eq. \2Q\ for the single impurity case. Difi^erent colours cor- 
respond to different localization regimes obtained from the 
solution of Eq. (|2ip . No real solutions of Eq. (|21[) correspond 
to the red region. If Eq. H21[) has real solutions, we can distin- 
guish three cases: no solutions are in the single particle spec- 
trum (yellow region), one solution is in the spectrum (green 
region), two solutions are in the spectrum (blue region). The 
four markers in the diagram correspond to the simulations 
presented in Fig. and to the values of indicated by ver- 
tical lines in Fig. [2] 



right shows a tendency to diverge due to the fact that the 
resonance lies just outside the single impurity spectrum. 

Let us recall that, in our single impurity analysis, we 
considered the case where one isolated dipole induces just 
a trimer of site energies {Ai,Ao,Ai} and we neglected 
beyond nearest neighbour contributions. In other words, 
we neglected the effect of long-range correlations. 

In the next section, we study in detail the role played 
by the dipolar tails that we neglected in this simplified 
calculation and we highlight the role played by long-range 
correlations. 



B. Long-range correlations 

In order to understand the role played by long-range 
correlations and place the dipolar case in a wider context, 
we investigate the localization properties of a disordered 
potential generated by an effective impurity-test particle 
interaction with tails decaying as u^^^ ^ l/l"-!^ where 
/3> 1. 

This is done by placing the impurities exactly as done 
in the dipolar case, keeping fixed the values of Aq and 
Al and choosing A„ = u^/ J = A/|n|'^ for n > 2. The 
case /3 = 3 recovers our physical model with dipolar in- 
teractions. Moreover, we shift and normalize the on-site 
energies in order to obtain the same average value (e„) 
and disorder strength W that we had in the dipolar case. 
Following this procedure, we can really analyze the ef- 
fect of long-range correlations keeping fixed the disorder 
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FIG. 5: (Color online) Localization length Cioc in units of 
the lattice spacing di, as a function of the energy in units 
of J for different type of long-range correlations identified by 
the asymptotic decay of the tails of the two-point correlation 
function C(£). The two different panels correspond to two 
different localization regimes induced by short-range correla- 
tions (see text for more details). 

strength W . In particular we considered values of P rang- 
ing from 1 up to 5 and we also considered the case of 
(3 = CO that corresponds to A„ = for |n| > 2. 

The potential generated with this procedure has Cp {£) 
which is unchanged and decays exponentially as previ- 
ously discussed. This is due to the fact that the impu- 
rities are placed exactly in the same way as before. The 
correlation function associated to the interaction poten- 
tial Cu {t) is instead modified and using Eq. one can 
show that it decays at large distances as C{i) ~ for 
/3 > 1 and as C[t) - log(^)/£ for (3 = 1. This asymptotic 
expressions determine also the shape of the tails of the 
reduced correlation function associated to the full poten- 
tial Ce{i). 

The effects played by long-range correlations are again 
studied by calculating numerically the localization length 
with the renormalization-decimation approach. In Fig. [S] 
we show the localization length Cioc calculated for dif- 
ferent values of /3. In particular we show a comparison 
between the two limiting cases oi (3 — oo and /3 = 1 
and the physical case under consideration, i.e. the dipo- 
lar case /3 — 3. We considered also other values of /3 
but we do not show the results here since they are not 
particularly instructive. They just show an intermediate 
behaviour between the two limiting cases reported here. 

The two panels in Fig. [5] correspond to two different 
set of parameters taken from Fig. [21 In the left panel we 
considered the case where the localization length is al- 
ways finite, while in the right panel the case where there 
are two resonances in the spectrum. Therefore the two 
curves for /3 — 3 are exactly the same curves which are 
shown in the upper row of Fig. [3] and those are compared 
with the case of complete absence of long-range corre- 
lations (/? = oo) and the case of very slowly decaying 
correlations (/? — 1). 



In Fig. [5l there are two main features that we would 
like to stress here. The first is the effect introduced by 
long-range correlations on the divergences that we have 
discussed in the previous section. We observe that they 
tend to be beveled. In fact in the 13 — oo case the lo- 
calization length takes values of the order of the system 
size N = 10'' signaling the presence of real metallic states 
in the system which extend over the full lattice. As the 
value of /3 is decreased the peaks in the localization length 
are still there but their height is decreased and they are 
shifted towards the band edges. This behaviour, that 
is reduced by lowering the impurity concentration C, is 
somehow expected, since the perfect resonance condition, 
obtained with the single impurity calculation presented 
in the previous section, was carried out neglecting the 
slowly decaying tails. When those tails are taken into 
account, they tend to restore interference effects in the 
scattering process and reintroduce localization in the sys- 
tem. However, let us stress that the effect of the short- 
range correlations remains clearly visible also in presence 
of long-range correlations. 

The second feature that we would like to highlight is 
the counterintuitive behaviour of Cioc introduced by long- 
range correlations at the center of the band. In fact, de- 
pending on the set of parameters under consideration, 
long-range correlations have an opposite effect on the lo- 
calization length. In the left panel, they introduce a de- 
crease of the localization length while in the right panel 
they tend to enhance it. This observation shows the 
highly nontrivial role played by long-range correlations 
in determining the localization properties of a disordered 
system and indicates a richer behaviour with respect to 
what has been observed so far in the literature [22| . 

These features that we extracted from the numerical 
simulations reported in Fig. [5] are also captured within 
the Born approximation. We do not report the curves for 

(2) 

^loc calculated in Born approximation, and just comment 
that the agreement between those curves and the exact 
numerical results is good, similar to what observed in 
Fig. [3) Finally we would like to comment that we do not 
find the presence of mobility edges induced by long-range 
correlations as suggested in j23l. |43|. 



V. CONCLUSIONS 

In this paper, we considered a set of dipolar impuri- 
ties pinned at random positions in a deep optical lat- 
tice which act as a disordered potential for an atom in 
a different internal state. An analysis of the statistical 
properties of the model showed that repulsive dipolar in- 
teractions between impurities introduce short-range cor- 
relations due to the fact that occupations of neighboring 
sites are forbidden. The localization properties of the 
model were calculated by means of a renormalization- 
decimation technique which allowed us to calculate prop- 
erties of very large systems and study the extended or 
localized nature of the states. We found that the pres- 
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cncc of short-range correlations can give rise to differ- 
ent regimes. In particiilar, as the parameters of the sys- 
tem are changed, wc observed regimes where one or more 
discrete sets of extended states appear in the spectrum. 
The occurrence of the different regimes can be predicted 
starting from an analytical expression obtained from the 
scattering of a single impurity. Long-range correlations 
were studied not only for the dipolar case but also for a 
more general two-point correlation function decaying as 
C{t) ~ l/^'^, where the case /3 = 3 corresponds to the 
dipolar case. We saw that long-range correlations tend 
to restore localization in the spectrum and lead to coun- 
terintuitive behaviours of the localization length. More 
precisely, depending on the regime under consideration, 
they can enhance or reduce localization at the center of 



the band. Our work sheds light on the interplay between 
the role of short-range and long-range correlations and 
can be a guide for experiments devoted to the study of 
Anderson localization with ultracold dipolar gases. 
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